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Abstract 

Weakly Interacting Massive Particles (WIMPs) are one of the leading candidates for 
Dark Matter. So far the usual procedure for constraining the WIMP-nucleon cross sections 
in direct Dark Matter detection experiments have been to fit the predicted event rate based 
on some model(s) of the Galactic halo and of WIMPs to experimental data. One has to 
assume whether the spin-independent (SI) or the spin-dependent (SD) WIMP-nucleus 
interaction dominates, and results of such data analyses are also expressed as functions 
of the as yet unknown WIMP mass. In this article, I introduce methods for extracting 
information on the WIMP-nucleon cross sections by considering a general combination 
of the SI and SD interactions. Neither prior knowledge about the local density and the 
velocity distribution of halo WIMPs nor about their mass is needed. Assuming that an 
exponential-like shape of the recoil spectrum is confirmed from experimental data, the 
required information are only the measured recoil energies (in low energy ranges) and the 
number of events in the first energy bin from two or more experiments. 



1 Introduction 



Astronomical observations and measurements indicate that more than 80% of all matter in the 
Universe is dark (i.e., interacts at most very weakly with electromagnetic radiation and ordinary 
matter) . The dominant component of this cosmological Dark Matter must be due to some yet to 
be discovered, non-baryonic particles. Weakly Interacting Massive Particles (WIMPs) x arising 
in several extensions of the Standard Model of electroweak interactions are one of the leading 
candidates for Dark Matter. WIMPs are stable particles with masses roughly between 10 GeV 
and a few TeV and interact with ordinary matter only weakly (for reviews, see Refs. [1, 2]). 

Currently, the most promising method to detect different WIMP candidates is the direct 
detection of the recoil energy deposited in a low-background underground detector by clastic 
scattering of ambient WIMPs off target nuclei [3, 4]. The basic expression for the differential 
event rate for elastic WIMP-nucleus scattering is given by [1]: 
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Here R is the direct detection event rate, i.e., the number of events per unit time and unit mass 
of detector material, Q is the energy deposited in the detector, po is the WIMP density near the 
Earth, ctq is the total cross section ignoring the form factor suppression and F{Q) is the elastic 
nuclear form factor, fi{v) is the one-dimensional velocity distribution function of the WIMPs 
impinging on the detector, v is the absolute value of the WIMP velocity in the laboratory frame. 
The reduced mass mr,N is defined by 

mr,N = ^ , 2) 

where is the WIMP mass and mN that of the target nucleus. Finally, Wmin is the minimal 
incoming velocity of incident WIMPs that can deposit the energy Q in the detector: 



^'min = ay (5 , (3) 

with the transformation constant 



and Vmax is the maximal WIMP velocity in the Earth's reference frame, which is related to the 
escape velocity from our Galaxy at the position of the Solar system, v^sc ^ 600 km/s. 



1.1 WIMP— nucleus cross section 

The total WIMP-nucleus cross section ao in Eq. (1) depends on the nature of WIMP couplings 
on nucleons. Generally, for non-relativistic WIMPs, one can distinguish spin-independent (SI) 
and spin-dependent (SD) couplings. 



1.1.1 Spin-independent couplings 



Through e.g., squark and Higgs exchanges with quarks, WIMPs could have a "scalar" interaction 
with nuclei^. The total cross section for the SI scalar interaction can be expressed as [1, 2] 

^o' = (^) <N [ZU + - Z)U] ' . (6) 

Here mr,N is the reduced mass defined in Eq. (2), Z is the atomic number of the target nucleus, 
i.e., the number of protons, A is the atomic mass number. A — Z is then the number of neutrons, 
/(p n) are the effective scalar couplings of WIMPs on protons p and on neutrons n, respectively. 
Here we have to sum over the couplings on each nucleon before squaring because the wavelength 
associated with the momentum transfer is comparable to or larger than the size of the nucleus, 
the so-called "coherence effect". 

In addition, for the hghtest supersymmetric neutrahno, and for all WIMPs which interact 
primarily through Higgs exchange, the scalar couplings are approximately the same on protons 
and on neutrons [5]: 

/n=:^/p. (7) 

The "pointlike" cross section ctq^ in Eq. (6) can thus be written as 

where mr,p is the reduced mass of the WIMP mass and the proton mass mp, and 

is the SI WIMP-nucleon cross section. The tiny mass difference between a proton and a neutron 
has been neglected. 

1.1.2 Spin-dependent couplings 

Through e.g., squark and Z boson exchanges with quarks, WIMPs could also couple to the spin 
of target nuclei, an "axial-vector" (spin-spin) interaction. The SD WIMP-nucleus cross section 
can be expressed as [1, 2]: 

= (f ) GWr,^ (^) [(^p)«P + (^n)aj' . (10) 

Here Gp is the Fermi constant, J is the total spin of the target nucleus, (5'(p,n)) arc the ex- 
pectation values of the proton and neutron group spins, and a(p,n) a-re the effective SD WIMP 
couplings on protons and on neutrons. 

For the SD WIMP-nucleus interaction, it is usually assumed that only unpaired nucleons 
contribute significantly to the total cross section, as the spins of the nucleons in a nucleus are 



^Besides of the scalar interaction, WIMPs could also have a "vector" interaction with nuclei [1, 2]: 



^vector ( ^ \ ™2 



2Z6p + {A- Z)bn , (5) 



where 6(p^„) are the effective vector couplings on protons and on neutrons, respectively. However, for Majorana 
WIMPs (x = x)) e.g., the lightest neutralino in supersymmetric models, there is no such vector interaction. 
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Table 1: List of the relevant spin values of the nuclei used for simulations presented in this 
paper. More details can be found in e.g., Refs. [1, 6, 7, 8]. 

systematically anti-aligned^. Under the "odd-group" assumption, the SD WIMP-nucleus cross 
section can be reduced to 

,SD /^32\ 2 _2 ('-^^^\lc \2| 



^o" = \- ) Gf<^ ) (%n))^|a(p,„)|^ . (11) 

Since for a proton or a neutron J = \ and (S'p) or i^S^^ — |, the SD WIMP cross section on 
protons or on neutrons can be given as 

'^xfp.n) = (^) G'FmJ(p_„)|a(p,n)|^ . (12) 

Moreover, once the upper and/or lower hmits on the WIMP-nucleon cross sections have been 
estimated by Eq. (12), it has been shown that, for a 'particular WIMP mass, one can use the 
following inequality to give constraints on the SD WIMP-nucleon couplings on the Op — On plane 
[6, 10, 8]: 

Qp ^ Qn \ ^ ^ ( ap ^ On \ 

Here 0-^^'^^^)'^''^'^''°^^'^^ are the upper/lower limits on the SD WIMP-proton/neutron cross sections, 
respectively, and the "±" sign in the parenthesis is the same as that of the (^n) / (Sp) ratio. So 
far the best constraint on the SD WlMP-proton coupling comes from the NAIAD [11], KIMS 
[12], SIMPLE [13], PICASSO [14], and COUPP [15] experiments: |ap| < 0.4 (for a WIMP mass 
of 50 GeV/c^) [13], whereas the best one on the SD WlMP-neutron coupling comes from the 
CDMS-II [16], XENONIO [17], and ZEPLIN-III [18] experiments: |a„| < 0.2 (for a WIMP mass 
of 50 GeV/c^) [13]. On the other hand, the relative strength of two couphngs for neutrahno 
WIMPs has been calculated as [19]^ 



0.55 < 



<0.8. (14) 



^However, more detailed nuclear spin structure calculations show that the even group of nucleons has some- 
times also a non-negligible spin (see Table 1 and e.g., data given in Refs. [1, 6, 7]). Hence, due to the neglect 
of the contribution from the even group of the target nucleons, the (exclusion limit of the) WIMP-nucleon cross 
sections could be overestimated. 

^However, Ellis et al. have shown that, in different theoretical scenarios, a„ could also be slightly greater than 
cip [20, 5]. 



Remind that the above conventional data analyses are independent of models of WIMP-nucleon 

couplings, but they do depend on the model of the Galactic halo through the use of the local 
WIMP density, po, and the velocity distribution of incident WIMPs, fi{v). Additionally, the 
results depend also strongly on the as yet unknown WIMP mass (see e.g., Refs. [10, 7, 8]). 



1.1.3 Comparison of the SI and SD interactions 

As discussed above, WIMPs could have both SI and SD interactions with target nuclei. Thus 
the WIMP-nucleus cross section gq in Eq. (1) should be a combination of the SI cross section iTq^ 
in Eq. (6) and the SD cross section in Eq. (10). However, due to the coherence effect with 
the entire nucleus shown in Eq. (8), the cross section for scalar interaction scales approximately 
as the square of the atomic mass number of the target nucleus. Hence, in most super symmetric 
models, the SI cross section for nuclei with /I > 30 dominates over the SD one [1, 2]. Nevertheless, 
as discussed in Refs. [21, 22, 23], in Universal Extra Dimension (UED) models, the SD WIMP 
interaction with nucleus is less suppressed and could be compatible or even larger than the SI 
one. 



1.2 Nuclear form factor 

1.2.1 For the spin-independent cross section 

For the SI cross section, there are some analytic forms for the elastic nuclear form factor. The 
simplest one is the exponential form factor, first introduced by Ahlen et al. [24] and Freese et 
al. [25]: 



Here Q is the recoil energy transferred from the incident WIMP to the target nucleus, 
1.5 



<5o = 



rriN-Ro 



is the nuclear coherence energy and 



Rq — 



0.3 + 0.91 



VGeV; 
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(15) 



(16) 



(17) 



is the radius of the nucleus. The exponential form factor implies a Gaussian form of the radial 
density profile of the nucleus. This Gaussian density profile is simple, but not very realistic. 
Engel has therefore suggested to use the following one [26], which derives from the nuclear 
density profile obtained by convolving a constant nuclear density with a gaussian one [27], and 
is similar to the numerical form factor derived from the Woods-Saxon nuclear density profile 
[1, 2], 

2 



fUQ) 



qRi 



Here ji (x) is a spherical Bessel function. 



q = \l2m^Q 
is the transferred 3-momentum, 



5s2 



(18) 



(19) 



(20) 



is the effective nuclear radius^ witfi^ 

Ra 1.2 A^/^ fm, (24) 

and 

s ~ 1 fm (25) 
is tlie nuclear skin thickness. 

1.2.2 For the spin— dependent cross section 

For the SD cross section, the form factor is different from nucleus to nucleus and no simple 
analytic form can provide a very good approximation. Generally, the form factor for the SD 
cross section can be expressed as [4, 1] 

Here the "spin structure" function S{q) depends generally on the SD WlMP-nucleon couplings: 
S{q) = alSooiq) + alSu{q) + aoaiS'oi(g) , (27) 

with the isoscalar and isovector coefficients: 

oo = Op + , ai = Op - On , (28) 

and 5*00 , -^n, and Sqi are the isoscalar, isovector and interference contributions to S{q), respec- 
tively. 

However, Klapdor-Kleingrothaus et al. have used the following form factor for the SD cross 
section [29], introduced by Lewin and Smith with the so-called "thin-shell" approximation [4]: 



ji(qRi) , for qR, < 2.55 or qR, > 4.5 , 

(29) 

const. ~ 0.047 , for 2.55 < qRi < 4.5 . 



1.2.3 Zero momentum transfer approximation 

For our simulations presented in this article, we will use the form factors given by Eqs. (18) 
and (29) for the SI and SD cross sections, respectively. However, it will be seen later that, 
since one would only have to estimate values of the form factors at the lowest energy ranges 
( < 20 keV for some currently running and projected experiments), we could practically use the 
"zero momentum transfer" approximation: 

F^{Q ~ 0) ~ 1 (30) 

in the methods introduced in this article. 



*In the literature, another often used analytic form for Ri has been given as [27, 4] 



i?i = y^^JT(ih2;^2^ , (21) 

where 

i?A ^ (1.23 A^/^ -0.6) fm, roc^0.52fm, s ~ 0.9 fm. (22) 

^For Ri given by Eq. (20) with s ~ 1 fm, another analytic form for R^ has also been given [28, 4]: 

Ra ^ (1.15 _^ 0.39) fm. (23) 



1.3 Motivation 



50 far the usual procedure for estimating the (exclusion hmits of the) WIMP-nucleon cross 
sections in direct Dark Matter detection experiments have been to fit the predicted event rate, 
dR/dQ in Eq. (1), based on some modcl(s) of the Galactic halo from cosmology and of WIMPs 
from particle physics to experimental data. Meanwhile, one has to assume whether the SI or 
the SD WIMP~nucleus interaction dominates. However, WIMPs should in general have both 
interactions with target nuclei. Moreover, as mentioned above, although in most models with 
neutrahno WIMPs as the best motivated candidate for Dark Matter, the theoretical predicted 

51 WIMP-nucleus cross section should be (much) larger than the SD one [5], Bertonc et al. have 
shown that another Dark Matter candidate, the lightest Kaluza-Klein particle (LKP) arising 
in the Universal Extra Dimension (UED) models, has a relatively larger SD cross section, or, 
equivalently, a larger cr^(^ to ratio [21]. Hence, for determining the nature of Dark Matter 
particles and distinguishing them between e.g., the lightest neutralino in supersymmetric models 
and the lightest Kaluza-Klein particles in models with Universal Extra Dimensions, estimates 
of both SI and SD cross sections, or, at least an estimate of the ratio between these two cross 
sections, in direct Dark Matter detection experiments is essential. 

On the other hand, as shown in our earlier work [30, 31] that one can determine the WIMP 
mass with direct Dark Matter detection experiments without a prior knowledge of the WIMP- 
nucleus cross section nor assumptions about the local density and the velocity distribution func- 
tion of halo WIMPs. It is therefore important to investigate methods for, conversely, extracting 
information on the WIMP-nucleon cross sections from experimental data without knowing the 
WIMP mass. 

The remainder of this article is organized as follows. In Sees. 2 and 3 I will show how 
to determine ratios of WIMP-nucleon couplings/cross sections once positive signals have been 
observed. Both the case that the SD WIMP interaction dominates (in Sec. 2) and that of a 
general combination of the SI and SD cross sections (in Sec. 3) will be considered. In Sec. 4 

1 will extend the data analysis procedure to the estimates of ratios between the SI WIMP 
scalar/vector couplings on protons and on neutrons. I conclude in Sec. 5. Some technical details 
for the data analysis will be given in an appendix. 

2 Only a dominant SD WIMP-nucleus cross section 

In this section I consider at first the case that the SD WIMP-nucleus interaction dominates over 
the SI one and derive the expression for determining the ratio between two SD WIMP-nucleon 
couplings. 

2.1 General expression 

By using a time-averaged recoil spectrum, and assuming that no directional information exists, 
the normalized one-dimensional velocity distribution function of halo WIMPs, /i(v), has been 
solved from Eq. (1) analytically [32] and, consequently, its generalized moments can be estimated 



by [32, 31]^ 
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r(gmin)/i^'(Qmin) + (n + l)/n(Qmin, Qr 



1 /'2 

2Qmin^(Qniin)/-^^(Qmin) "I" -^o(Qmin7 Qmax) 



(31) 



Here ^(Q) = oi.\/Q^ Q(min,max) ^le the experimental minimal and maximal cut-off energies of the 
data set, respectively, 

\ ^ / expt, (5=(9min 

is an estimated value of the measured recoil spectrum {dR/dQ)expt {before normalized by an 
experimental exposure S) at Q = Qmm, and In{Qmin, Qmax) can be estimated through the sum: 



-^n(Qmin) Qmax) ^ ' 



Aftot Q{n-l)/2 



11 Qv 



(33) 

and A^tot is the 



where the sum runs over all events in the data set that satisfy Qa € [Q-d 
number of such events. 

Now, since the integral on the right-hand side of Eq. (1) is just the minus-first moment of the 
velocity distribution function, {v^), which can be estimated by Eq. (31), by setting Q — Qmin 
and using the definition (4) of a, one can obtain straightforwardly that 



2 



(34) 



F (Qmin) 

Then, in order to avoid the uncertainty of po (of a factor of ~ 2 [1]), one can combine two 
experimental data sets with different target nuclei, X and Y, to eliminate po in Eq. (34) and 
thus obtain the following expression for the ratio between the WIMP cross section on nuclei X 
and Y: 

^ = ^,/^f^V (35) 

where mr,(x,y) are the reduced masses of the WIMP mass and the masses of target nucleus, 
"^(x,y), and I have defined 



1 



2Qmin, Xl^x{Qmm,x) 
Fx{Qmm,x) 



o,x 



(36) 



and similar for TZ^^^y F(^x,y) iQ) form factors of the nucleus X and Y, r(x,Y) {Qmin,(x,Y) ) 

refer to the counting rates for the target X and Y at the respective lowest recoil energies in- 
cluded in the analysis, and £{x.y) are the experimental exposures with the target X and Y. 
The emphasize here is that Eq. (35) can be used once positive signals are observed in two (or 
more) experiments; information on the local WIMP density po and on the velocity distribution 
function of halo WIMPs, fi{v), are not necessary'''. 



^Here we have implicitly assumed that Qmax is so large that terms involving — 2(5max r((5max)/-F^((5max) 
are negligible. Due to sizable contributions from large recoil energies [32], this is not necessarily true, especially 
for some not very-high Qmax in the experimental reality, and/or heavy detector targets, and/or heavy WIMPs. 
Nevertheless, we will show in this and the next sections that, since we use only n = —1, 1, and 2, Eq. (31) 
can still be used for the determinations of the ratios between different WIMP-nucleon couplings/cross sections. 
Moreover, considering the large statistical uncertainties due to (very) few events in the highest energy ranges, 
this should practically be a good approximation. 

''Later we will see that nor information on the WIMP mass is necessary. 



Substituting the expression (10) for ctq^ into Eq. (35) and using the definition (4) of a for 
both target nuclei, one can solve the ratio between two SD WIMP-nucleon couphngs analytically 
as [33] s 



i^f" {s^)x±{s^)Y{nj,n,xinj,n,Y) 

for n 7^ 0. Here I have used the relation [31]: 

Ox _ 'R'n,Y 
OtY T^n,X ' 

and defined 
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Jx + ly7^n,x 

with TZcr,x defined in Eq. (36) and 



(38) 



(39) 
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■n,X 



2QLty Vx(Qn.in,x)/Fl(Qmin,x) + (^ + l)In,X 



2Qmin,x'^x(Qmin,x)/-^x(<5min,x) + h,X 



l/n 



(40) 



T^j,n,Y and 7^„,y can be defined analogously^. Note that Eq. (37) is independent of the WIMP 
mass and can be used for estimating a^/a-p with measured recoil energies directly. 

Because the couplings in Eq. (10) are squared, we have two solutions for On/op here; if exact 
"theory" values for TZj,n,{x,Y) are taken, these solutions coincide for 
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for nj^n,x = , 



for 7^J,„,y = , 



(41) 



which depends only on properties of two used target nuclei (sec Table 1). Moreover, it can be 
found from Eq. (37) that one of these two solutions has a pole at the middle of two intersections, 
which depends simply on the signs of {S,^)x and (5*11)^: since TZj^n,x and 7lj^n,Y are always 
positive, if both (S'n)x and (S'n)y are positive or negative, the "— (minus)" solution (an/ap)^^„ 
will diverge and the "+ (plus)" solution {a^/ttp)^^ will be the "inner" solution; in contrast, 
if the signs of (5'n)x and {S^)y are opposite, the "— (minus)" solution (on/op)^'^^ will be the 
"inner" solution (sec Figs. 1). 

SD 

By using the standard Gaussian error propagation, the statistical uncertainty on (an/ap)_,_^ 
can be expressed as 




{Sp)Y{Sa)x — {Sp)x{Sa)Y 



1 7^ 



J,n,X 



{Sn)x ± ('S'n)y (7^ J,n,x/'^J,n,y) 



2 TZjnV , 



^Note that, although the constraints on two SD WIMP-nucleon couplings have conventionally been shown in 
the flp — On plane, considering the theoretical expected value given in Eq. (14), we use always the a^/a^ ratio in 
our work. 

^Hereafter, without special remark all notations defined for the target X can be defined analogously for the 
target Y and occasionally for the target Z. 



1 f dnn,x \ _ j_ f dn,,x \ 

T^n,X V ) T^a,X \ dCi^x ) _ 

' 1 ( &llr,,x \ 1 ( dlla,x W , . 

\ 1/2 

+ (X ^ F) I . (42) 

Here a short-hand notation for the six quantities on which the estimate of (cin/ap)±J^n depends 
has been introduced [31]: 

C\,X = In,X , C2,X = Io,X , Cz,X = rx{Qmm,x) ] (43) 

and similarly for the Q^y. Estimators for cov(ci, Cj) and explicit expressions for the derivatives 
of 7ln,x and TZ^^x with respect to q x will be given in the appendix. Note that 7?.o-,(x,y) are 
actually independent of Ci^(x,y) = In,{x,Y)^ n ^ 0. 

In Figs. 1 I show the numerical results for a target combination of ^^Ge and ^^Cl with 5,000 
experiments based on the Monte Carlo simulation^°. The theoretical predicted recoil spectrum 
for the shifted Maxwellian velocity distribution [1, 2, 32] with a Sun's orbital velocity in the 
Galactic frame vq = 220 km/s, an Earth's velocity in the Galactic frame Vc = 1.05 vq/^ and a 
maximal cut-off velocity of the velocity distribution function v^ax = 700 km/s, as well as the 
nuclear form factor given in Eq. (29) have been used. The experimental minimal and maximal 
cut-off energies have been set as Qmin = 5 keV and Qmax = 100 keV for both targets. Each 
experiment contains an expected number of 50 total events; the actual event number is Poisson- 
distributed around this expectation value. The input WIMP mass has been set as 100 GeV. 

As discussed above, since (<S'n)73Ge and (5*11)3701 have the same sign, the "+" solution shown 
in the left frame of Figs. 1 is the inner solution for the range of interest < a^/ap < 1, while 
the "— " solution shown in the right frame diverges between — (5'p)73Ge/('S'n)73Ge — —0.08 and 
— (5'p)37ci/(5'n)37ci = 1.16. Note here that, for practical use of analyzing real data, one might 
however not be able to make a choice from the "+" and "— " estimates given by Eq. (37), 
especially if they are close to the coincidences, e.g., around 1.16 or —0.08 here. For example, for 
a true Cn/ap = 1.1, one will get (an/ap)^ = 1.1 and (an/ap)^ = 1.25, the same results as for 
the case with a true a-^/ap — 1.25. 




2.2 Reducing statistical uncertainty on (an/ap)^^„ 

For estimating the statistical uncertainty on (an/flp)|°„ by Eq. (42), one needs to estimate 
contributions from the counting rate at the threshold energy, r(Qmin), from given in Eq. (33), 
and from the covariance between r((5min) and /„. From Eqs. (A9), (AlO) and (A13) in the 

appendix, one can find a way to reduce these statistical uncertainties by estimating the counting 
rate, instead of at the experimental minimal cut-off energy, at the shifted point Qs^i (from the 
central point of the first bin, Qi): 

Qs,i = Qi + 7- In 



^"Note that, rather than the mean values, in this article we give always the median values of the reconstructed 
results from the simulated experiments. 

^^The time dependence of the Earth's velocity in the Galactic frame [1, 2] has been ignored. 



sinh(/i;i6i/2) 



hbi/2 



(44) 




Figure 1: The reconstructed a^ja^ ratios estimated by Eq. (37) and the lower and upper bounds 
of their la statistical uncertainties estimated by Eq. (42) with n = — 1 (dashed blue), 1 (solid 
red), and 2 (dash-dotted cyan) as functions of the input an/ctp ratio. Here I show the "+ (— )" 
solutions in the left (right) frames separately. The theoretical predicted recoil spectrum for the 
shifted Maxwellian velocity distribution with = 220 km/s, Vc = 1.05 fo, and v^aax = 700 km/s 
as well as the nuclear form factor for the SD cross section given in Eq. (29) have been used. 
^■^Ge and ^^Cl have been chosen as two target nuclei. Each experiment contains 50 total events 
on average in the energy range between 5 and 100 keV. The input WIMP mass has been set as 
100 GeV. See the text for further details. 



where ki is the logarithmic slope of the reconstructed recoil spectrum in the first Q— bin and bi 
is the bin width. Then, according to Eq. (A9), the measured recoil spectrum at Q = Qs i can be 
estimated by 

6i ' 




r{Qs,i) = [ ;^ ] 

expt, 1, Q=Qs,i 

with the statistical uncertainty given as 



(45) 



a 



[riQs,i)) 



- (n) = ^ 



(46) 



where A^^i is the event number in the first bin. 

In Figs. 2 I show the reconstructed a^/a^ ratios and the lower and upper bounds of their 
lo" statistical uncertainties with n = —1 (dashed blue), 1 (solid red), and 2 (dash-dotted 
cyan) estimated by Eq. (37) with the counting rates at the shifted points of the first Q— bin, 
'^{x,y)(Qs,i,(x,y)) = ^\x,Y),i as functions of the input a^/a^ ratio"^^. It can be seen that the sta- 
tistical uncertainties on [a^/a^)^^ estimated with different n (namely with different moments 



of the WIMP velocity distribution) with r(^x,Y){Qs,i,{x,Y) 
almost equal. Therefore, since 

1/2 
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) are clearly reduced and, interestingly. 
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^^Labeled hereafter with an "sh" in the subscript. 



'^Ge + ^^Cl, Q„|„ > 5 keV, Q^^ < 100 keV, 2 x 50 events, = 100 GeV "Ge + ^'Cl, Q„in > 5 keV, Q^,,^ < 100 keV, 2 x 50 events, m,, = 100 GeV 




Figure 2: As in Figs. 1, except that we estimate (an/ap)|^„ with the counting rates at the shifted 
points of the first Q-bin, r(x,y)(Qs,i,(x,y)) = ?^(x,y),i- 

one would practically only need events in the lowest energy ranges (~ 20 events between 5 
and 15 keV in our simulations) for estimating a^/a^. Consequently, one has to estimate the 
values of form factors only aX Q = Qs,i, and the zero momentum transfer approximation 
F'^{Q ~ 0)) ~ 1 can be used. In fact, our simulation shows that a relatively higher thresh- 
old energy (Qmin ~ 10 keV and Qs^i ~ 14 keV) should not affect the reconstruction of an/ap 
significantly, especially for the first approximation with pretty few events and thus a large sta- 
tistical uncertainty. 

On the other hand, as mentioned above, the expression (37) for estimating the ratio be- 
tween two SD WIMP-nucleon couplings is independent of the WIMP mass. In Figs. 3, I show 
the reconstructed an/ap ratio and the lower and upper bounds of their la statistical uncer- 
tainties as functions of the input WIMP mass for a fixed input a^/ap = 0.7. We estimate 
with r(x,y)(Qmin,(x,y)) and r(^x,Y)(Qs,i,(x,Y)) ^^^^ right frames, respectively. It can 

be seen that, firstly, except the statistical uncertainty estimated with r(x,y)(<5min,(x,y)) and 
n = —1 (the dashed blue curves labeled as (an/ap)rec, so, -i in the left frame), for WIMP masses 
> 50 GeV, the reconstructed an/ap ratio as well as the statistical uncertainty are (almost) 
independent of the WIMP mass; however, if WIMPs are (very) light {m^ < 25 GeV), a^/ap will 
be (strongly) underestimated, due to the non-zero threshold energies^^. Secondly, the statistical 
uncertainties on a^/ap estimated with r(x,y)(<5min,(x,y)) and r(x,y)(<5s,i,(x,y)) are only 10% or 
even 7% combined with an ~ 1.5% systematic deviation. 

As a comparison, I show the combinations of the "+" and "— " solutions with n = 1 shown 
in Figs. 1 and 2 together in the left frame of Figs. 4. In the right frame, I compare also the 
results with n = 1 shown in Figs. 3. The ~ 30% (from 10% to 7%) reduction of the statistical 
uncertainty by estimating with r(x,y)(<5s,i,{x,y)) for > 100 GeV can be seen obviously. 

-"^^Remind that, as discussed in Ref. [34] for the method for estimating the SI WIMP-nucleon couphng, this 
kind of underestimate (or overestimate shown later in this article) can be alleviated (corrected) once we can 
decrease the threshold energies (to be negligible); see also Ref. [35] for simulations with negligible experimental 
threshold energies. 




Figure 3: The reconstructed a^ja^ ratios estimated by Eq. (37) and the lower and upper bounds 
of their la statistical uncertainties with n = — 1 (dashed blue), 1 (solid red), and 2 (dash-dotted 
cyan) as functions of the input WIMP mass m^. Here we estimate with r(x,y)(Qmin,(x,y)) (left) 
and r(x,y)(<5s,i,(x,y)) (right). The input a^/a^ ratio has been set as 0.7. The other parameters 
and notations are as in Figs. 1 and 2. 

Furthermore, considering the low natural abundances of ^'^Ge and '^^Cl (see Table 1), in 
Figs. 5 we simulate with another combination of target nuclei: ^^F and ^^''I. As discussed in the 
previous subsection and shown in Figs. 1, 2, and 4, the inner solutions of {a^/a-^)^^ have a much 
smaller statistical uncertainties and the range of these inner solutions depends on the —{Sp)/ {S^) 
values of our target nuclei. Hence, one benefit of using the combination of ^^F and ^^''I is that 
one can estimate (an/ap)|°„ in a much wilder range of interest: |an/ap| < 4. Consequently, for 
the practical use of analyzing real data, one has therefore not to worry about making the choice 
from the "+" and "— " estimates, which is discussed at the end of the previous subsection; since 
{Sn)wF and (S'n)i27i have different signs, we can just take the "— " solution in Eq. (37). 

However, Figs. 5 show us also some drawbacks of the use of the ^^F + ^^''I combination. 
For WIMP masses > 50 GeV, (an/ap)|°„ estimated with r(x,y)(<5min,(x,y)) (dashed blue) 
are ~ 15% - 30% overestimated; whereas those estimated with r(x,y)(Qs,i,(x,y)) (solid red) are 
even worse: ~ 36% for ~ 1 TeV. Moreover, the statistical uncertainties shown here become 
also much larger (of a factor of ~ 3 — 5) than those shown in Figs. 4. This enlargement of the 
statistical uncertainties is mainly caused by the larger value of the prefactor of cr((an/ap)|^„) in 

Eq. (42). According to Table 1, the values of (5'p)y(S'n)x — ('S'p)x('S'n)y are ~ 0.023 for the 
Ge + CI combination, but ~ 0.067 for F + I. Meanwhile, as shown in both Figs. 4 and 5, the 
statistical uncertainties are at the largest in the middle of two coincidence points and reduce 
as the input a^/a^ approaches to one of these two points. Since we set the input On/ap = 0.7 
for simulations with different WIMP masses, comparing to the relative difference between 0.7 
and the middle point of 1.16 and —0.08, i.e., 0.54, the relative difference between 0.7 and the 
middle point of 4.05 and —4.12, i.e., —0.035, is slightly smaller. This causes also a larger 
statistical uncertainty for the use of the F + I combination. In contrast, for light WIMP masses 
[nT'x ~ 50 GeV), the two estimates with ^^F and ^^^I shown in the right frame of Figs. 5 are 
much more underestimated than results shown in the right frame of Figs. 4. 




Figure 4: Left: the combined results of the "+" and "— " solutions estimated with 
i^{x,Y)iQrain,ix,Y)) (dashed blue) shown in Figs. 1 and with r(^x,Y)iQs,i,ix,Y)) (sohd red) in Figs. 2. 
Right: a comparison of the results shown in Figs. 3 estimated with r(x,y)(Qmin,(x,y)) (dashed 
blue) and with with r(x,y)(<5s,i,{x,Y)) (solid red). Only the results estimated with n = 1 are 
shown here. 




Figure 5: As in Figs. 4, except that ^^F and ^^^I have been chosen as two target nuclei. Note 
that the scales used here are different from those in Figs. 4. 



3 Combination of the SI and SD cross sections 



In this section I consider the general combination of the SI and SD WIMP-nucleus cross sections. 



3.1 General expression 

At first, by combining Eqs. (8), (10), and (12), we can find 
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For the general combination of the SI and SD WIMP-nucleus cross sections, the expression (1) 
for the differential event rate should be modified to 
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where I have used Eq. (8) again. Then one can find straightforwardly that the integral above 
can be estimated by Eq. (31) with the following replacement: 

(51) 



Hence, for this general case, Eq. (34) becomes to 
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(52) 



(53) 



Now by combining two targets X and Y and using the definition (4) of a, the relation (38) 
between ax/cnY with n = —1, as well as the expression (40) for 7^n,(x,y), one can obtain that-*^^ 



-^SI,x(Qmin,x) + 



^Xp/^Xp) ^P'^-^SD,x(Qmin,x) TZ 



-^SI,y(Qii 



,y) + (</^xp) Cp.y^lD,y(Qmin,y) ^-.^^ ' 



(54) 



^^This equation can be obtaincsd by simply assuming that the integral over fi{v) on the right hand side of 
Eq. (50) estimated in two experiments (approximately) agree and can thus be cancelled by each other. This 
assumption can practically always hold, even though the experimental minimal and maximal cut-off energies 
in these two experiments should be matched by requiring [31] that a.x yJQ(m\i\,ma-x).x = <^y -\/Q(min,max).y ■ 
since, as the expressions (56) and (57) show, only the estimated values of ri^x,Y){Qmm,{x,Y)) are important 
for the data analysis. Note that, however, once one applies similarly this simple cancellation for the case of 
a dominant SD WIMP cross section discussed in the previous section, only the expression (37) for [a^/ a^)^^ 
with n = —1, namely with TZj^-i.(x,Y) given in Eq. (47), can be obtained. This is because that, by using 
this cancellation, a cx ^m-^/mj-^-^ and 2(5^j^r((5min)/-F^('5min) + -^o on the right-hand side of Eq. (34) will 
be eliminated before one obtains this equation. Then one cannot use the relation (38) to convert ax/ ay to 
T^n,Y /Ti-n,x and therefore to obtain the expression (37) with different values of n; except with n = — 1, since 

x''' x{Qvain,x) / F'^[Qjaij^^x) + Io,x appears in the numerator oi TZ-i^x (see Eq. (40), not in the denominator 
as for the cases with n = 1, 2, • • •) and can thus be cancelled out anyway. 



where I have assumed m(x,Y) oc >1(x,y) and defined 

^-^= Exml ■ ^^^^ 

From Eq. (54), the ratio of the SD WIMP-proton cross section to the SI one can be solved 
analytically as [33] 

"^XP _ -^SI,y('5min,y)('^m,x/'^m,y) " ^^SI,xiQmm,x) ^^^^ 



m,y J 



where Cp^(x,y) have been defined in Eq. (49). Similarly, the ratio of the SD WIMP-neutron cross 
section to the SI one can be given analogously as [33]^^ 

_ -^SI,y(Qmin,y)(7^m,x/7^m,y) - -^SI.X (Qmin.x) ^^^^ 

C^XP Cn,xF^B,x{Qmm,x) " Cn,y-FsD,y (Qmin.y) ('7^TO,x/'7^m,y ) ' 

with the definition 

' (^p)(«p/«n) + (^n) l\ ^5g^ 

The emphasize here is that one can use expressions (56) and (57) to estimate c''x(p,n)/'^xp '^'^^^out 
a prior knowledge of the WIMP mass m^. Moreover, since C(p,n),(x,y) depend only on the nature 
of the detector materials, o"^(^n)/^xp practically only functions of 'R-m,{x,Y), i-e., the counting 
rates at the experimental minimall cut-off energies, which can be estimated by using events in 
the lowest available energy ranges. 




3.2 Using (a„/ap)SD gq. (37) 

Since Cp and C„ defined in Eqs. (49) and (58) are functions of an/^p, once the fln/op ratio has been 
estimated (from e.g., some other direct detection experiments by Eq. (37) under the assumption 
of a dominant SD WIMP-nucleus interaction), cr^p /a^p can then be estimated by Eq. (56) with 
the following statistical uncertainty 
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Explicit derivatives of c^p /c^p with respect to Cp^(x,y) and TZjn,{x,Y) will be given in the appendix. 
Note that Eq. (37) can be used only when the SD WIMP-nucleus interaction really dominates 
over the SI one. We will see later that, if the SD interaction does no^ dominate, the a-a/a^ ratio 
should not be estimated by Eq. (37) any more. 



^^Herc I assumed that ct^Jj ~ cr^p by Eq. (7). 
Hereafter I consider only the case with protons. But aU formulae given in this section can be applied 
straightforwardly to the case with neutrons by replacing p — >■ n and Cp C^. 



3.3 Solving a^/ap with a third nucleus 



Nevertheless, for the general combination of the SI and SD WIMP-nucleus cross sections, the 
On/op ratio can in fact be solved analytically by introducing a third nucleus with only an SI 
sensitivity: 
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I.e., 



Cp,z = . (62) 
Then, according to Eq. (56), we have 

z{Qmin,z){T^m,X /T^m,z) " -^fl.X (Qmin.x) -^SI,z(Qmin,z) ("T^m.y /"T^m.z) " -^SI,y (Qmin,y ) 
Cp,xFs^x{Qnim,x) Cp,y Fl^y (Qmin.y ) 

Using Cp defined in Eq. (49), the On/op ratio can be solved analytically as [33] 
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Here I have defined 
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Note that, firstly, (an/ap)^^^^ and Cp,(x,y) given in Eqs. (63), (64a), and (64b) are functions 
of only r(^x,Y,z){Qmin,{x,Y,z)), which can be estimated with events in the lowest energy ranges. 
Secondly, while the decision of the inner solution of (an/ap)±^n depends on the signs of (5'n)x and 
(S'n)y, the decision with {(hjap)^^^^ depends not only on the signs of Sn/p,x = {Sn)x/ {Sp)x and 
Sn/p,y = (>S'n)y/(>S'p)y, but also on the order of the two targets. For the Ge + CI combination, 
since Sn/pjsQ^ = 12.6 > Sn/p 37ci = —0.86, one should use the upper expression in the second line 
of Eq. (63), and since Sn/p/sge and Sn/p,3rci have the opposite signs, the "— (minus)" solution 
of this expression (or the "+ (plus)" solution of the expression in the first line) is the inner 
solution. In contrast, since Sn/p,i9F = —0.247 < Sn/p,i27i = 0.243 and since Sn/p,i9F and Sn/p,i27i 
have the opposite signs, the "— (minus)" solution of the lower expression in the second line of 



Eq. (63) (or the "— (minus)" solution of the expression in the first hne) is then the inner solution 
for the F + I combination. 

Finally, from the expression (63), the statistical uncertainty on {an/ ap)^'^^^ can be given by 
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And the statistical uncertainty on the ratio between two WIMP-proton cross sections in Eq. (56) 
can be expressed as (c.f., Eq. (59)) 
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with d (o-Jp/o-Jp) /5(a„/ap) given in Eq. (60) and 
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for i — X, y, Z. Explicit derivatives of (on/op)!^''"^^ and Cp^(x,Y) ^^^^ given in the appendix. 

In Figs. 6 I show the reconstructed a^/dp ratios estimated by Eqs. (37) (dashed blue) and 
(63) (solid red) and the lower and upper bounds of their la statistical uncertainties estimated 
by Eqs. (42) and (66) as functions of the input a^/ap ratio^^. For the SI cross section the nuclear 
form factor given in Eq. (18) has been used. The SI WIMP-proton cross section has been set as 
10~^ pb (left) and 10"^° pb (right), respectively, whereas the SD WIMP-proton couphng Op has 
been set as 0.1.^^ Besides '''^Gc and ^^Cl, ^^Si has been chosen as the third target for estimating 
Cp,{x,Y) ^'i^- (64a) and (64b). 

In the left frame, it can be seen obviously that a^/ap estimated by Eq. (37) (dashed blue) 
under the assumption of a dominant SD WIMP-nucleus interaction has two discontinuities 
around (an/ap)in = 1-16 and —0.08 and the reconstructed a^/ap ratio is systematically over- 
/underestimated. In contrast, On/ap determined by Eq. (63) (solid red) shows a more smooth 
estimate, although the reconstructed ratio is a bit underestimated with a relatively larger sta- 
tistical uncertainty for input a^/ap ratios around 1.16. However, once we set the input SI 
WIMP-proton cross section two orders of magnitude lower and thus the SD WIMP-nucleus 
cross section really dominates (the right frame), the On/op ratios estimated by two methods 
show a clear compatibility. 

In Figs. 7 the first two targets with both SI and SD sensitivities have been replaced again by 
^^F and ^^^I. In contrast to Figs. 6, a^/ap estimated by Eqs. (37) (dashed blue) and (63) (solid 



^^Note that all results shown in this subsection are only reconstructed with T(^x.y.z){Q s,i.(x,y.z)) = '"(x,y,z).i- 
^^Remind that the current exclusion limit on the SI WIMP-nucleon cross section is < 5 x 10"^ pb for WIMP 
masses of ~ 30 - 100 GeV from the XENONIO [36], CDMS-II [37], XENONIOO [38], and EDELWEISS-II [39] 
experiments, whereas the limits on the SD WIMP couplings on protons and on neutrons are |ap| < 0.4 and 
[an I ^ 0.2 (for a WIMP mass of 50 GeV/c-^) [13], respectively. On the other hand, the theoretically predicted 



values for is |ap| 



< 



0.1 [5]. 




Figure 6: The reconstructed a^/a^ ratios estimated by Eqs. (37) (dashed blue) and (63) (sohd 
red) and the lower and upper bounds of their la statistical uncertainties estimated by Eqs. (42) 
and (66) as functions of the input a-a/a-p ratio. Besides ''^Ge and ^^Cl, ^^Si has been chosen as the 
third target for estimating Cp,(x,y) by Eqs. (64a) and (64b). For the SI cross section the nuclear 
form factor given in Eq. (18) has been used. The SI WIMP-proton cross section has been set 
as 10~^ pb (left) and 10"^° pb (right), respectively, whereas the SD WIMP-proton coupling Op 
has been set as 0.1. The other parameters are as in Figs. 4. Note here that the scales of the 
(an/ap)rec^axes in two frames are different. See the text for further details. 

red) shown here are overestimated, especially the ratio reconstructed under the assumption of a 
dominant SD interaction. Nevertheless, the an/op ratio estimated by Eqs. (63) (solid red) in both 
Figs. 6 and 7 show that the ratio between two SD WIMP-nucleos couplings could in principle 
be estimated correctly with an ~ 20 — 40% statistical uncertainty without prior information on 
the WIMP mass nor on the SI WIMP-nucleon cross section. The (in) compatibility between 
the reconstructed a^/ap ratios under different assumptions and/or with different combinations 
of target nuclei could also allow us to check whether the SD WIMP-nucleus interaction really 
dominates or not. 

Similar to the right frames of Figs. 4 and 5, Figs. 8 show the reconstructed ctn/c^p ratios 
estimated by Eqs. (37) (dashed blue) and (63) (solid red) and the lower and upper bounds of 
their la statistical uncertainties estimated by Eqs. (42) and (66) as functions of the input WIMP 
mass m^. The over-/underestimated a^^/ap ratios with different combinations of target nuclei 
can be seen obviously here. For input WIMP masses < 50 GeV, all estimates are as usual 
(strongly) underestimated. Nevertheless, for WIMP masses > 50 GeV, the reconstructed 
la statistical uncertainty intervals estimated by Eqs. (63) and (66) (solid red) could basically 
cover the input (true) value pretty well. 

3.4 Choosing nuclei with C^^y = and {Sp)x ^ {Sn)x — 

In the expression (56) for the ratio of two WIMP-proton cross sections, there are four sources 
contributing statistical uncertainties, i.e., Cp^(x,Y) T^m,{x,Y)- order to reduce the statistical 
uncertainty on the estimate of cTyp/Cyp, one can choose at first a nucleus with only an SI 
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Figure 7: As in Figs. 6, except that ^^F, ^^^I, and ^'^Si have been chosen as target nuclei. The 
scales used here are also different. 



'^Ge + ^'Cl (+ '^''Si), Q^in > 5 keV, Q^^^ < 100 keV, 2(3) x 50 events, ( 
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Figure 8: The reconstructed an/op ratios estimated by Eqs. (37) (dashed blue) and (63) (solid 
red) and the lower and upper bounds of their la statistical uncertainties estimated by Eqs. (42) 
and (66) as functions of the input WIMP mass m^. The input an/op ratio has been set as 0.7, 
the other parameters are as in Figs. 6 and 7. Left: ^^Ge, ^''Cl, and ^^Si have been chosen as the 
three target nuclei. Right: ^^F, ^^''I, and ^^Si have been chosen. Note here that the scales of the 
(an/ap)rec^axes in two frames are different. 
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Figure 9: The reconstructed o'^p/c^p (left) and cr^n/'^xp (^ig^^) the lower and upper bounds 
of their la statistical uncertainties as functions of the input a^/ap, respectively. The dashed 
blue curves indicate the values estimated by Eq. (56) with a^/ap estimated by Eq. (63) {not by 
Eq. (37)); whereas the solid red curves indicate the values estimated by Eq. (71). ^^Ge has been 
chosen as the second target with only an SI sensitivity and combined with ^^Na (for a^/a^p) 
and ^^^Xe (for a^/a^p) for using Eq. (71). Parameters are as in the left frame of Figs. 7, the 
input SI WIMP-proton cross section has been set as 10~^ pb. Note that, since we fix cx^p and Op, 
the theoretical curve of a^/a^p is a constant, whereas the curve of a^/a^p oc a parabola. 
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Figure 10: The reconstructed o'^p/'^xp (1^^^) '^xn/'^xp (^^§1^^) lower and upper bounds 

of their la statistical uncertainties as functions of the input WIMP mass m^, respectively. The 
input SI WIMP-proton cross section has been set as 10~^ pb and the input a^/ap = 0.7, the 
other parameters and notations are as in Figs. 9. 



sensitivity as the second target: 



(^p)y = (^n)y = 0, 



(69) 



I.e., 

Cp,y = 0. 

The expression in Eq. (56) can thus be reduced to [33] 
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Then we choose a nucleus with a (much) larger proton group spin as the first target: 

(^p)x»(^n)x^O, 

in order to eliminate the an/op dependence of Cp,x given in Eq. (49)^^: 
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and the statistical uncertainty given in Eq. (67) can be reduced to 
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(75) 



In Figs. 9 I show the reconstructed c^p/c^p (left) and cr^J^/cr^p (right) as functions of the 



XP 



Xn ' 



XP 



input On/op, respectively. The dashed blue curves indicate the values estimated by Eq. (56) with 
fln/flp estimated by Eq. (63) {not by Eq. (37)); whereas the solid red curves indicate the values 
estimated by Eq. (71). Since, as shown in the left frame of Figs. 7, fln/op can be estimated pretty 
well by Eq. (63) with the target combination F and I in the range of interest —1 < |an/ap| < 2, 
cr^p / shown in the left frame here can be reconstructed with an ~ 40% statistical uncertainty 
by using the combination of ^^F + ^^^I + ^*Si. On the other hand, the right frames of Figs. 9 
and 10 show also that o'^n/o'^p could still be estimated well with -^^-^Xe and '''^Ge by Eq. (71), 
although the statistical uncertainty is now larger (~ 70%). 



4 Estimating ratios of the SI WIMP— nucleon couplings 

So far I have used the theoretical prediction (7) that the SI scalar WIMP couphng on protons 
is approximately equal to the coupling on neutrons. For the sake of completeness, I consider in 
this section briefly the case that WIMPs have different SI scalar or vector couplings on protons 
and on neutrons [40]. For WIMPs having only the scalar interaction with nuclei, the expression 
(6) for (Tq^ can be rewritten as 
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n 2 



(76) 
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One can also choose {S^x ^ {Sp)x — and C^^x given in Eq. (58) becomes 
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Thus one can obtain the following replacements: 



J+1 
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and 



Z 
I' 



Z 



Substituting Eq. (77) into Eq. (39), we can get 
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(77) 
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(79) 



where 7^o-,x and TZn,x are given in Eqs. (36) and (40).^° Then the ratio between the scalar 
WIMP coupling on protons and on neutrons can be estimated analogously to Eq. (37) as 



7n 



Zx ± ZY\/'R-fT,x/T^n,X\/T^n,Y/'^cr,Y 
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with the following statistical uncertainty: 
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(81) 



Note that, firstly, since A — Z > for all nuclei, the inner solution of /n//p given in Eq. (80) with 
a much smaller statistical uncertainty is always the "+" solution. Secondly, the two coincident 

points of the "+" and "— " soulutions decided by —Zxl{Ax — Zx) and —ZyI{Ay — Zy) arc 
however always negative. While, for lighter nuclei, e.g. ^^Si and ^^F, the values of —Zj {A — Z) 
are ~ —1; for heavier nuclei, e.g. ^^^I or ^^^Xe, these values are ~ —0.7. This means that, 
unfortunately, for confirming the /n//p ratio with the theoretical predicted value of ~ 1, we can 
only use the "outer (— )" solutions given in Eq. (80) with much larger statistical uncertainties 
and data sets with piles of events should therefore be required. 

On the other hand, assuming that WIMPs have only the vector interaction with nuclei, 
according to the expression (5) for o"o*^'^*°'', we can write down the expression for the relative 
strength of two "vector" couplings directly as 



X\ 2Zx ± 2ZY^n,,x/nn,xyJn^,Y/n,,Y 
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(82) 



^•^Remind that the form factor F'^{Q) here must be chosen for the SI cross section. 



with the following statistical uncertainty: 
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Note that the factor "2" appearing in the denominator of the prefacor in Eq. (81) has been 
cancelled here. 



5 Summary and conclusions 

In this paper, I presented methods for determining ratios between different WIMP-nucleon cou- 
plings/cross sections from elastic WIMP-nucleus scattering experiments. All methods presented 
here are independent of the model of halo WIMPs as well as (practically) of the as yet unknown 
WIMP mass. Assuming that an exponential-like shape of the recoil spectrum is confirmed 
from experimental data, the required information arc only the mcasm^cd recoil energies and the 
number of events in the first energy bin from at least two direct detection experiments with 
different detector materials having spin sensitivities contributed from protons and/or from neu- 
trons. Even better, our simulations show that, for estimating the relative strengths of different 
WIMP-nucleon couplings, one would only need events in the lowest available energy ranges. 

In order to avoid the uncertainty on the local WIMP density po, our analyses are based 
on combining two (or more) experiments using different target nuclei. By assuming, as the 
first step, that the SD WIMP-nucleus interaction dominates over the SI one, the expression 
for determining the ratio between two SD WIMP-nucleon couplings, ttn/ap, has been rederived 
[33]. Then our simulations with different combinations of target nuclei show that, in order to 
obtain an unambiguous result with much smaller statistical uncertainty in the range of interest: 
lon/flpl ^ 2, nuclei with sensitivities on both protons and neutrons should be more suitable than 
nuclei being sensitive (almost) only on protons or on neutrons. 

More generally, I considered also the combination of the SI and SD WIMP-nucleus cross 
sections. By using three different targets, two of them have non-zero group spins from pro- 
tons and/or from neutrons, the second expression for determining the ratio between two SD 
WIMP-nucleon couplings can be rederived [33]. Although its statistical uncertainty depends on 
the relative strength between the SD and SI WIMP-nucleus interactions, the (in) compatibility 
between the a^/ap ratio reconstructed under different assumptions and/or with different combi- 
nations of target nuclei could allow us to check whether the SD WIMP-nucleus interaction really 
dominates. Moreover, by using two or three different nuclei, one or two of them have non-zero 
group spins from protons and/or from neutrons, one can in principle also determine the ratios 
of the WIMP-proton/neutron cross sections to the SI ones, o"^p/o"^p and o"^n/o-^p, directly. 

Our simulations presented here are based on several simplified assumptions. Firstly, the 



sample to be analyzed contains only signal events, i.e., is free of backgrounds^' Secondly, 
all experimental systematic uncertainties as well as the uncertainty on the measurement of the 
recoil energy have been ignored. The energy resolution of most currently running and projected 
detectors is so good that its uncertainty can be neglected compared to the statistical uncertainty 
with (very) few events in the foreseeable future. 

In summary, I demonstrated in this paper the use of our new methods for extracting infor- 
mation on WIMP-nucleon couplings/cross sections, which are independent of models of WIMPs 
from particle physics as well as of models of the Galactic halo from cosmology. By combining 
with information on the estimation of the SI WIMP-nucleon coupling [44, 34], one could in 
principle estimate the absolute values of the spin-dependent couplings/cross sections. These 
information could help us not only to give constraints on different models of particle physics in 
the parameter space, but also to understand the nature of halo Dark Matter particles as well as 
to distinguish them between candidates predicted in different scenarios [21, 22, 23, 5]. 
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A Lists of needed formulae 

Here I hst all formulae needed for our model-independent data analyses described in this article. 
Detailed derivations and discussions can be found in Refs. [32, 31]. 



A.l Estimating r(Qmin) and /n(Qmin, Qmax) 

First, consider experimental data described by 

Qn-'-f <Qn,^<Qn+'-f, i = 1 , 2, ■ ■ ■ , U = 1 , 2 , ■ ■ ■ , B . (Al) 

Here the total energy range between Qmin and Qmax has been divided into B bins with central 
points Qn and widths 6„. In each bin, A„ events will be recorded. Since the recoil spectrum 
dR/dQ is expected to be approximately exponential, the following ansatz for the measureo? recoil 
spectrum {before normalized by the experimental exposure £) in the nth bin has been introduced 
[32]: 

fdB\ = rne^-{Q-Q.,r.) _ (A2) 

) expt, n J expt, Q~Q„ 

Here r„ is the standard estimator for {dR/dQ)e^pt at Q — Qn. 

rn = ^ , (A3) 



^^For background discrimination techniques and status in currently running and projected direct detection 
experiments see e.g., [41, 42, 43, 37], 

^^For detailed simulations and discussions about effects of residue background events on the determinations of 
ratios between different WIMP couplings/cross sections see [35]. 



kn is the logarithmic slope of the recoil spectrum in the nth Q— bin, which can be computed 
numerically from the average value of the measured recoil energies in this bin: 
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The error on the logarithmic slope kn can be estimated from Eq. (A4) directly as 
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Qs,n in the ansatz (A2) is the shifted point at which the leading systematic error due to the 
ansatz is minimal 
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(A8) 



Note that Qs,n differs from the central point of the nth bin, Qn- Prom the ansatz (A2), the 
counting rate at Q — Qmin can be calculated by 

r(Qn.in) = rie'=^(«-'"-«-i), (A9) 

and its statistical error can be expressed as 



(j'(r(Qn.in)) = r\Q, 
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Finally, since all are determined from the same data, they are correlated with 

A'tot Q{n+m-2)/2 

cov(4,/J = E^^, 
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(A12) 



where the sum runs over all events with recoil energy between Qmin and Qmax- And the correlation 
between the errors on r((5min); which is calculated entirely from the events in the first bin, and 
on In is given by 
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note that the sums li here only count in the first bin, which ends at Q — Qmin + 

On the other hand, with a functional form of the recoil spectrum (e.g., fitted to experimental 

data), {dR/ dQ) cxpt, one can use the following integral forms to replace the summations given 
above. Firstly, the average Q— value in the nth bin defined in Eq. (A5) can be calculated by 
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For IniQmin, Qmax) givcu in Eq. (33), we have 

/•Qmax Q("-1)/2 /'dR\ 

and similarly for the covariance matrix for 7„ in Eq. (A12), 

r-Qma. g{n+m-2)/2 /^^\ 



(A14) 



(A15) 



cov(/„, I„ 
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Remind that {dR/dQ)cxpt is the measured recoil spectrum before normalized by the exposure. 
Finally, /^(Qmin, Qmin + ^i) needed in Eq. (A13) can be calculated by 
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Note that, firstly, r(Qinin) and /« (Qmin i Qmin + ^i) should be estimated by Eqs. (A9) and (A17) 
with ri, ki and Qs^i estimated by Eqs. (A3), (A4), and (A8) in order to use the other formulae for 
estimating the (correlations between the) statistical errors without any modification. Secondly, 
r{Qmin) and /n(Qmin, Qmax) estimated from a scattering spectrum fitted to experimental data 
are usually not model-independent any more. Moreover, for the use of Eqs. (33), (A12), (A15), 
(A16), and (A17) the elastic nuclear form factor F^{Q) should be understood to be chosen for 
the SI and SD WIMP-nucleon cross section correspondingly. 



A. 2 Derivatives of 'T^n,x and IZ^^x 

First, from Eq. (40) one can find explicit expressions for the derivatives of TZn,x with respect to 
clx are: 
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explicit expressions for the derivatives of TZn,Y with respect to Cj,y can be given analogously. Note 

that, firstly, factors TZn,{x,Y) appear in all these expressions, which can practically be cancelled 



by the prefactors in the bracket in Eq. (42). Secondly, all the /| 



0,{X,Y) 



and /, 



n,{X,Y) 



should be 



understood to be computed according to Eq. (33) or (A15) with integration limits Qmm and 
Qmax specific for that target. 

Similarly, expressions for the derivatives of TZ^^x can be computed from Eq. (36) as 
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and similarly for the derivatives of TZa,Y- Remind that factors 7?-(T,(x,y) appearing here can also 
be cancelled by the prefactors in the bracket in Eq. (42). 



A. 3 Derivatives of crt^/<7^\ 

From the expression (56) for estimating cr^p/c^p, its derivatives with respect to Cp^(^x,Y) '^^^ 
given as 
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Meanwhile, the derivatives of c^p/'^xp ^^^^ respect to TZrn,{x,Y) are 
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On the other hand, from expression (49) for Cp one can find that 
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A.4 Derivatives of (an/ap)^^+^° 

At first, from the first and second lines of the expression (63), we have, 
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Then, from the definition (64a) of Cp^x, one can get directly 
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Similarly, from the definition (64b) of Cp,y, we have 
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